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Abstract 

An effective intra- and inter-ladder charge-spin hamiltonian for the quarter- 
filled ladder compound a'-NaX^Os has been derived by using the standard 
canonical transformation method. In the derivation, it is clear that a finite 
inter-site Coulomb repulsion is needed to get a meaningful result otherwise 
the perturbation becomes ill-defined. Various limiting cases depending on the 
values of the model parameters have been analyzed in detail and the effective 
exchange couplings are estimated. We find that the effective intra-ladder 
exchange may become ferromagnetic for the case of zig-zag charge ordering 
in a purely electronic model. 

We estimate the magnitude of the effective inter-rung Coulomb repulsion 
in a ladder and find it to be about one-order of magnitude too small in order 
to stabilize charge-ordering. 



PACS: 71.10.Fd Lattice Fermion models (Hubbard model etc.), 75.30.Et Ex- 
change and super-exchange interactions, 64.60.-i General studies of phase 
transitions. 
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I. INTRODUCTION 



Low- dimensional quantum spin systems have received considerable attention from both 
theoretical as well as experimental point of view due to their unconventional physical prop- 
erties. a'-NaV205, which was believed to be a low- dimensional inorganic spin-Peierls (SP) 
compoundB has recently been under intense investigation. ct'-NaX^Os is an insulator and its 
magnetic susceptibility data fits very well to the one-dimensional Heisenberg chain model 
yielding an exchange interaction J=440 and 560 K for temperatures below and above the 
spin-Peierls transition temperature Tsp (Tsp~ 34 K ) respectivelyBI. For T< Tsp, an 
isotropic drop in the susceptibility corresponding to a singlet-triplet gap of Asp=85 K has 
been observed. 

Recent X-ray structure data analysis!'! at room temperature disfavours the previously re- 
ported non-centrosymmetric structure! (Cj v — P2imn space group) where V +4 spin- 1/2 ions 
form a one- dimensional Heisenberg chain, running along the crystallographic b-direction, 
separated by chains of V +5 spin- zero ions. But the evidence for the centrosymmetric point 
group! (-D2I — Pmmn) leads to only one type of V-site with a formal valence +4.5 in this 
compound. The V-sites then form a quarter-filled ladder, running along the b-axis with the 
rungs along the crystallographic a-axis. In the quarter-filled scenario, the electron spins are 
not localized at V-ions rather distributed over a V-O-V molecule which has found support 
by NMPiH as well as Raman measurements!. 

The nature of the state below Tsp is presently under intense investigation. Isobe and 
UedJ originally proposed a usual spin-Peierls scenario but the detection of two inequivalent 
V-sites in nmr! indicates a more complicated scenario and the possibility of charge ordering. 
Several types of charge ordering, including 'in-line'! and 'zig-zag ordering has been 
proposed, but only the zig-zag type of ordering has been found to be in agreement with 
neutron scatter ing00 and anomalous X-ray scattering^!. 

Recent determinations of the low-temperature crystal structure found the space group 
Fmm20§ and proposed the existence of three inequivalent V-ions below This 
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scenario was investigated by DMRG (density-matrix renormalization group) and a cluster- 
operator theory0 and a strong disagreement with neutron scattering dataEil was found. 
Ohama et al. recentely observed^ that the aparent contradiction between cyrstallography 
(three inequivalent V-ions bewlow Tsp) and NMR (two inequivalent V-ions bewlow Tsp) 
could be resolved when one considers possible subgroups of the originally proposed space 
group Fmm2lli. 

The crystal structure of a'-NaX^Os at T>T sp is orthorhombic (a=11.318 A, b=3.611 A, 
c=4.797A) and consists of double chains of edge-sharing distorted tetragonal VO5 pyramids 
running along the orthorhombic b-axis. These double chains are linked together via common 
corners of the pyramids and form layers. These are stacked along c-direction with no-direct 
V-O-V links. The Na atoms are located in between these layers. For the orbitals of the 
d-electrons at V-sites, those with d xy symmetry are suggested to be the relevant ones above 
and below Tgp§. Due to the special orbital structure, the hopping amplitudes t a and i& are 
much larger than the inter-ladder hopping t ab , t a and tb being the hopping amplitudes along 
the rung and the ladder direction respectively (t a ~ 0.38 eV, tb ~ 0.17eV, t ab ~ 0.012 eV). 
Since a'-NaX^Os is an insulator, it has been assumed that the on-site Coulomb repulsion 
U is sufficiently large in comparision to the hopping amplitudes (U ~ 2.8 eV from DFT 
calculationB). Moreover, one has to introduce the inter-site Coulomb repulsions, V a , H, 
V a b to obtain the required charge ordering. In fact, it has been shown in a Hartree-Fock 
calculation^ that the condition U > V a , Vb, V a b > t a , t b , t a b must be fulfilled in order to achieve 
a complete charge ordering. We consider this and other limits in the present paper. 

In the present work, we take into account the charge dynamics to obtain an effective 
low energy hamiltonian for a'-NaX^Os. One starts from a pure electronic hamiltonian, 
which includes electron hopping in and between the ladders as well as the on-site and inter- 
site Coulomb interactions. The on-site Coulomb interaction U is taken to be the largest 
parameter in our calculation. Since a'-NaX^C^ is an insulator and we work at quarter-filling, 
one can project on a subspace of states which contains one electron per rung. Therefore, 
it is convenient to use an Ising pseudo-spin variable r 2 = ±1/2 corresponding to a rung 
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with an electron on the right/left site of the rung. This is in the same spirit of Kugel and 
Khomskii's treatment of the orbital degeneracy problem in Jahn- Teller systemslll. The spin 
and the pseudo-spin operators can be written as, 

s z = W< R l R ° + L] Mi s + = r\r 1 + l\l 1 , s- = r\r^ + l\l^ (1) 

r z = \YMR* - L l L ^ T+ = £^> = j: 4^, (2) 

Z a a a 

which for example yields, R^L^ = t?{\ + fif), iZ-.T-kl = rfSf, Lj 4 L i>T = (± - t?)%~, 
etc., where = r x + ir y , S' ± = S x + iS' i ' and R\ a (Li ia ) are the creation operator of an 
electron with spin a on the right (left) site of the i-th rung of the ladder. In (|l|) and (0) we 
have suppressed the site-indices. 

II. INTRA-LADDER EXCHANGE 

Let us start with an electronic hamiltonian for the quarter-filled ladder (see Fig. |T]), 
which can be written as, H — Hq + H' Q + Hj, with 

H = t a Y^{R\,a L i,<r + h.C.) + U Y^( n i,m n hRA + 1H,L,Hni,L,i) + V a ^2 n i,R,v n i,L,cr> , (3) 
i,cr i i,a,a' 

Ho = Vb X! ( n i,L,ani + l t L,a' + ^i,_R,cr^i+l,_R,cr' ) + V b " ^ ( n i,L,ani + l t Ry + n it R^rii + i t L,a') 

1,(7,(7' %,(7,u' 

+V a b n mi R i(T n n> L,o> , (4) 

<m,n> ,<7,cr' 

Hi = tb y^XRi.gRi+ha + L i>a Li+i )a + h.C.) + t a b ^ ^m,a^n,a + L^Rm^) , (5) 

i,o <m,n>,<7 

where t a , U and V a are the hopping integral, on-site and the inter-site Coulomb repulsion in 
a rung respectively, tb, V b ' and V b " are the hopping integral and the Coulomb interaction in 
between rungs in a ladder, whereas t ab and V ab are the inter-ladder hopping and the Coulomb 
interaction respectively. n^R^n^i^) is the electron density operator with spin a in the right 
(left) site of i-th rung and < m,n > denotes the pair of rungs m and n on adjacent ladders. 
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We estimate the parameters of the inter-site Coulomb repulsion using a screened Coulomb 
repulsion, V = e 2 /(ed), where e is the dielectric constant and d the distance between the 
respective Vanadium atoms. The distance in between two V-ions along the rung and the 
leg of the ladder (a- and b-directions) are 3.502 A and 3.611 A. The dielectric constant is 
e = 11 from microwave and far infrared measurements^. One obtains, V a = 0.3738 eV and 
V b ' = 0.3625 eV. The diagonal V-V distance in the b-direction in the ladder is 5.030 A, which 
implies, V h n = 0.2643 eV. It will be clear from the later discussion that the effective inter- 
rung Coulomb repulsion V& is given by the difference between V b ' and V b ", i.e. V b = V b — V b , 
which comes out to be small (V& = 0.1023 eV) compared to V a . Next, V a b = 0.4305 eV, as the 
V-V inter-ladder distance is 3.0401 A. Note that V ab is slightly higher than V a and nearly 
four times higher than Vb. 

In order to develop a perturbation expansion, we start by considering the case of a single 
two-leg quarter-filled ladder. The one-electron eigenstates of Hq for a single rung consist of 
bonding and anti-bonding wave functions, which we denote as a^|0 >= ~ L')\0 > and 

s^O >= 77|C^ + ^)|0 >, with eigenenergies —t a and +tb respectively. Now let us consider 
the coupling of the rungs along the legs described by the first term in the hamiltonian Hj. 
In order to obtain a coupling between the pseudo-spin and the spin variables, we use here 
the standard canonical transformation method^, which is given by, 

H eff = e !S i7e- iS , (6) 

where the operator S is determined from the condition 

Hj + i[S,H Q ] = 0, (7) 

which turns out to be 

$ = T^, =r\\n><n\Hi\n' ><n'\. (8) 

n,ri ~ ^n) 

Thus, the effective hamiltonian can be written as, 
H ef f = H + H' n - - V f — 1 - + - 1 - I In >< nlHAn' >< n'lHAn" >< n"\, 

" ° K^n»W-En E' n -E») 1 ' 71 ' 11 

(9) 
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where the initial and the final states \n" > and \n > are the two-rung states, i.e., all possible 
combinations of the bonding and the anti-bonding states between the nearest neighbour 
rungs. In the present case, there are sixteen possible such states which are the following: 

4.4-' 1 >> 4- a L'i° >> a U4>'i° >> a U a L' 1 >. ( 10 ) 

with a, a' =],[■ The six intermediate states \n' > which are the two-particle excited states 
in a rung, have to be antisymmetric under the exchange of both spin and pseudo-spin 
coordinates, in accordance with the Pauli principle. Thus, we have two sectors for the 
excited states depending on the total and the z-component of spin as well as the pseudo- 
spin quantum numbers which are labeled as, \S,S z ;r,r z >. Hence, the states involved 
are, |0,0;1,1 >= i^XjO >, |0, 0; 1, -1 >= 4t4J° >> I°>°;1,0 >= ^lAl ~ 
R[ l Ll A )\0> and |1,1;0,0>=4 T 4 T |0>, |1, -1; 0, >= ^LjjO >, |1, 0; 0, >= 
-^(RI^lI^ + R\^L\^) |0 >. The eigenenergies of the excited states in the large U limit 
are V a for the spin-triplet states \1,S Z ;0,0 > (S z = —1,0 + 1). For the spin-singlets, the 
eigenenergies are U' for ^(|0, 0;1,1 > +|0, 0;1, — 1 >) (symmetric), U for ^(|0, 0;1,1 > 
— 10, 0; 1, — 1 >) (antisymmetric) and V' a for |0, 0; 1, >, with U' w U + and V' a w 

K — u_ a v ■ After some lengthy but straightforward algebra and in the case of large but finite 
U, the total effective hamiltonian can be written as, 

H eff = H + H>+ HtfJ a + Hl n f f, (11) 

where B. l ^j a is the effective intra-ladder hamiltonian which one can express as, H l ^j a = 
Hfjf +Heff, with an< ^ Heff De i n g the contribution due to the intermediate spin-triplet 
and spin-singlet states respectively. H l JjJ r is the effective inter-ladder hamiltonian which 
can be derived in a similar way and will be discussed in the next section. In terms of the 
pseudo-spin and spin variables, the unperturbed hamiltonian H and H' can be expressed 
as, 

H = 2t a J2r*, (12) 

i 
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h> = 2(K - vi') E (7 + \ z t? + i) + n v h » + v ab E (t - « + ^V^) ' ( 13 ) 

i V ^ ' <m,n> V ^ & / 

where N is the number of rungs. In a similar way, the effective hamiltonian H^Jt can be 
written as, 



It is obvious from the above expression that Hf ff is independent of the Coulomb correlation 
energy U. This is due to the fact that while deriving this effective hamiltonian we have 
used the eigenenergies for the excited states which happen to be V a for this case. Since the 
effective hamiltonian H^h is obtained due to the contribution from the same intermediate 
states \n' > and different initial and final states |n > and \n" >, it can be written as, H^X = 
H^tf + #e// + ^eff (which are the contributions due to the antisymmetric, symmetric and 
1 0, 0; 1, > intermediate states), with 

Hiff =-^EQ- 2rf + n • f m ) (1 - S, ■ S l+1 ) , (15) 



2fl Efl- !t ^ 1 - 2 ^i+^-^ (t-^-^i), (16) 



C/' + 2t a V V4 2 * i+1 ' ,_r 7 V4 



^--^e(^^-^ 1+ -^)Q-^-^ 

2t * E (7 - - ^1 + ■ ^1) (7 - £ ■ • (17) 



v;' + 2t a y V4 2 1 4+1 1 lT v V4 

It should be noted here that one gets a non-zero contribution to H^jf and even if 

U = 00 which will be discussed below. Moreover, it is obvious from the expression for H^L 
(see Eq. (|TJ|)) that a finite V a is indeed needed to get a meaningful result otherwise the 
perturbation becomes ill-defined for V a = 0. 
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Limiting Cases and Discussion 

Case I: t a = 0: This limit implies U' = U and V' a = V a and thus, H^ff and -f^// can 
be combined together to yield, 

— ( ' 2, --fE(j + ^ 1 ) (i 



+ s '» , (18) 



whereas the effective hamiltonian H^f) gets reduced to, 

H ( eff = E (l + * ■ n + i - 2*7^) Q - Si • S i+1 ) . (19) 

Since H^jf neither depends on t a nor on U, it doesn't get affected in the above limiting case 
and the same is true for the other cases considered below. The effective hamiltonian derived 
by Thalmeier and FuldJ corresponds to Eq. (p~8|) . 

Case II: t a = 0, U — oo: In this limit, which also corresponds to the limit V a ^> 2t a , 
the contribution to the effective hamiltonian from -f^// and n^fl vanish and thus, the total 



intra-ladder effective hamiltonian becomes the sum of nflf an d Heff ( see Eq. (HID an< ^ 



(P~9|) ) . This is what is exactly obtained by Mostovoy and Khomskii0 but with a different 
interpretation. 

(si) (s2) 

Case III: U = oo, V a <C 2t a : In this case, the effective hamiltonian H e tj and H^J 
vanish but H^f) reduces to, 

H [ :Pf = f a E (*f + ^i) Q - ^ ■ SWi) , (20) 

so that the total effective intra-ladder hamiltonian becomes the sum of Eq. ( |T3| ) and (^U|) . 
Moreover, since V a <C 2t a , the major contribution will be from Eq. ([14]). 

Case IV: Disordered Phase: In the disordered phase, where the electrons are in the 
bonding states (see Fig. [| (a)), we can get an estimate of the effective exchange coupling 
in the effective hamiltonian by taking the averages over its charge (pseudo-spin) part. One 
can write down the effective exchange hamiltonian (disregarding the constant factors) as, 
^eff h = Jeff h ^iSi ' ^i+i- m the present case, we have, < rf >= —1/2 and < r\ >= 



=< rf > whereas < rfrf +l >= 1/4 and < rfrf +1 >= =< rfrf +1 >. Thus, the effective 
exchange coupling due to the hamiltonian H^h and H^lf vanish, but that of and 
i?e// become, jj,+£ t anc ^ v+2t wn ich yields Jfff 1 = 2t% ( v ,+2t + v+2t )• ^ * s c l ear that 
the exchange coupling here is antiferromagnetic (> 0). Using the parameters mentioned in 
the present work, Jfff 1 is estimated to be 0.08 eV. The expression for Jfff 1 is exactly the 
same (for the case V a = 0) as obtained by Horsch and Mack0. 

Case V: Complete Charge Ordered Phase: We can get an estimate of the effective 
exchange couplings in the completely charge ordered (zig-zag) phase (see Fig. |^ (b, c)), 
following the same procedure as that of the disordered phase. Here, one has, < rf >= 
=< rf > and < rf >= 1/2, < rf +1 >= -1/2 whereas < rfrf +1 >= =< rfrf +1 > and 
< rfrf +1 >= —1/4. Hence, the effective exchange coupling due to H^ff an d -^i// van ish 
but that of H®j and H^jj survive, which ultimately leads to Jtjf 1 = —2tl[y — 2(v-2t ) — 
However, using the parameter values, it is calculated to be -0.17 eV. In the case 



where U — oo and for < {V a /2t a ) < 1, Jfff 1 becomes ferromagnetic (< 0) whereas for 
(V a /2t a ) > 1, it is antiferromagnetic. The variation of Jftf 1 with respect to the parameter 



2(V a '+2t a )J 

/~v/"i fivin rr \w II 11/ t J-r I I , 

'eff 

r l /~vn rvr / 

eff 

(V a /2t a ) is shown in Fig. [|. It is clear from the figure that there exists a minimum in the 
ferromagnetic region for (V a /2t a ) = 1/V3, where J™f min = -2.26 (t 2 h /t a ) = -0.2 eV, which 
is quite large. 

III. INTER-LADDER EXCHANGE 

Next, let us consider the hopping between the two nearest neighbour ladders, which is 
described by the second term of the hamiltonian Hj (see Eq. (|5])), i.e., 

H} nter = t a b X/ (-^m,o-^n,o- + L n(T R mt(T ). (21) 

<m,n>,(T 

The effective inter-ladder coupling between the charge and the spin degrees of freedom is 
derived in the same way as has been done for the single ladder case. Since the two-particle 
excited states in this case are exactly the same as what has been done earlier, the effective 
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inter-ladder hamiltonian can be written as sum of two parts, i.e., H$f = if*Jf (t) + #^. er(s) , 
where the superscript and V have the same meaning discussed in the previous section. 
The effective hamiltonian H™^ r ^ is derived to be, 



2(K + 24) <S> 2 « + ^-^)(l + ^-^) 



E (l-2« + f m .f n )(^ + 5 ? m .5 ? „), (22) 

V a <m,n> V ^ / \4 / 

1 TT inter(s) 1 ... TT inter(s) TT inter(sl) , TT inter(s2) , TT inter(s3) 1 

whereas -n e // can be written as -n e // = n e jj + -" e // + e// > with 

— t 2 

jjinter(sl) _ L ab ST^ 

^ ~ 2{U - 2t a ) £ n> 

1 X + X Z _ 2 v , 1 

to 1 n _i_ to n O-7- 2 o- 2 t-^ o- 2 i t- 2 t- 21 i J z? \ / _ C 

^ ~r 2 "r 2 ZT m T n T m T n ' 'm'n ~T~ T m ' T n J I ^ "m ' 

2(C/ + 2£ a ) <m ^ n> 

1 x + x z _ z x / J 

to 1 ra 1 'to n r> 2 2 , Z z 2 a; , =? ^ | Q C 

^ 2 1 2 ZT TO T n T m T n T m T n T m ' T n J It ' 

4 < s>G +2i f s - ? -- f -) (!-*••«•)■ (23) 

The expression for H l ^ r ^ is exactly the same as that of H™^ 7 ^ except that U here is 
replaced by U' . On the other hand, _£/~*™^ er(s3 ) j s obtained as, 

rjinter(s3) _ tjb \- f 1 , T to + T n q y y , -? -Wl o .0 

**eff ~ 2(K'-2U < ^>U + 2 ^ T «+r»r B j^ 6 m 6 r 



+2 / 1 r x 1 T x \ / 1 _ _ 



2(K' + 2t a ) <^> V4 2 ™ » V V4 

e. ^ n . . . . . .\ /1 



-^77 E ( i -2« + r m -r ri )( i -5 ? m -5 ? n). (24) 



Since the derivation of the effective inter-ladder hamiltonian proceeds in the same way as 
that of the intra-ladder case, the limiting cases will follow the same way as has been done 
before. Moreover, here also, one needs a finite V a (but V a ^ 2t a ) in deriving the effective 
hamiltonian, otherwise the perturbation becomes ill-defined for V a = 0. Furthermore, one 
gets a non-zero contribution to the effective inter-ladder hamiltonian even if U = 00. 
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Limiting Cases and Discussion 



Case I: t a = 0: In this limit U' = U and V' a = V a follow naturally. Thus, H l JJ t j T ^ 
reduces to, 



Hi t;r{t) = x + ^ , + ^ _ & x (25) 



y a <m,n> V4: 7 
n . .1 i TT inter(sl) TT inter(s2) ^ TT inter(s3) i 

bimilariy, H e ft , H e ff an( ^ "eff become, 

Trinter{sl) Tjinter{s2) _ 4t afc ST^ , T m ~ T n z z\ q q \ /r> K \ 

U eff +tL eff - Jj- 2^ 2 m U J \4 _ ' J ' 

^e,(,3) = _2^ ^ (I +7^(1/4 (27) 

The expression Eq. (26) is the same one as has been obtained by Thalmeier and Fuldei. 

Case II: t a — 0, U — oo: In this case, the contributions from H l jff Sl ^ and H % ^ T ^ 
vanish. The contributions from //™** er w anc i ^nter(sz) can ^ e combined together to yield, 

rrinter Tjinter{t) rjinter{s) 2t afe ^ / 1 \ , , 

which becomes independent of the spin degrees of freedom. 

Case III: U = oo, V a <C 2t a : In this case, the contributions from H % ^ r ^ and H % ^ T ^ 
vanish. Combining H™^ rl ^ and H t e j t j V ^ s3 \ the spin-dependence drops out and one obtains, 

j.2 j.2 /i x 

Tj-inter _ Tj-inter(t) n inter (s3) _ ^ab ^ , an ''ab I Q-r 2 t x -U -r t I 

^e// ~ n eff +tL eff ~7T~ 2-^i \ m ' T n ) ^7 L 7 ZT m T n + T m T n I ■ 

™a <m,n> v a <jn,n> VH: 7 

(29) 

Case IV: Disordered Phase: Following the procedure already mentioned for the intra- 
ladder case, we can get an estimate of the effective exchange coupling between the nearest 
neighbour ladders (see Fig. || (a)) by taking the averages over the charge part in the ef- 
fective hamiltonian. Here also, one can write down the effective hamiltonian (disregard- 

jinter,exch __ / jinter,exch(t) jinte', 
l eff - We// + J eff 

'eff 
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,i , , r , \ rrinter^exch / T inter,exch(t) rinter,exch(s)\ ^ a a 

mg the constant factors) as, rl eff = {J e ff + J e ff )/^< m ,n> • o n 

jinter,exc J2 <mn> S m • S n . In the present case, one has, < >=< t% >= —1/2 and 



T m T n >• Thus, the effective exchange coupling which is due to /f^t , H^ r ^ sl \ H z e nt ' 

anrl £T inter ( s3 ) ;„ o-ivpn fl^ jinter,exch _ ~t\ b r i i i mi 

anu n ff is given db, J eff — „ L rv „, , (tja-*?.*- \ (TT>±->.t..\ <V4-?j-,\\- ex ~ 



< >=< >= =< >=< > whereas < r^r^ >= 1/4 and < r^T# >= =< 

T-mter(sl) rj-inte', 

t e// id givcu era, ^ e// - — L(v a +2t a ) ~ (!7+2t a ) ~~ (!7'+2t a ) ~ (V a '+2t a )] 

change coupling here is antiferromagnetic and is estimated to be (j™ t jr r > exch j '^2j = 0.39. For 
V a = 0, it gives rise to the same expression as obtained by Horsch and MackH. 

Case V: Complete Charge Ordered Phase: Here, we have four different completely charge 
ordered (zig-zag) phase depending on the state through which we compute the averages over 
the charge part of the effective hamiltonian. In all these cases, one has, < r^j >=< t% >= 
=< >=< > and < r^r^' >= =< t^t^ >. In addition to this, one has, 



to the state .R^L^ JO > (Fig. |3| (b)). The effective exchange coupling for ff 1 ™ 1 '"' 



(i) < >= 1/2, < >= —1/2 and < Tm r n >= — 1/4, where the averages are due 

jinte'i 
l eff 

and ff™y r ( s3 ) vanish but that of H™ff r and jj l ^ r ^ s2 > are finite which gives rise to, 
jintpexch _ ^ [_L_- _|_ + jj + 2 (u>-2t a ) + 2(u'+2t a ) + M • The exchange coupling 
here is antiferromagnetic and is estimated to be {J^f t f r,exch / t 2 ab ) = 1.42. For t a = 0, it 
gives rise to usual super-exchange, i.e. j l JJ t jr^ exch = 4tl b /U. 



(ii) < t;^ >= —1/2, < >= 1/2, and < t^t* >= —1/4, where the averages are taken 

jinte'i 
'eff 

r,exch 

'eff 



over the state L^J^JQ > (Fig. | (c)). In this case, j^ e ^« = = " M) 
which ultimately corresponds to the case j™^ r > exch = q. 



(iii) < >=< >= 1/2 and < t^t* >= 1/4, where the averages here are due to the state 
Rln R^plO > (Fig. || (d)). The effective exchange couplings for H™ff Sl ' and H % ^ T<yS2 ^ 

• u rni jinter,exch • 1 . j r .i . .-i j , TT inter(t) j rrinterfsS) 

vanish, thus, J e ^ is obtained from the contribution due to ti e ^ w and rl e ^ 

which is jinter,exch _ -ij b r l 1 i J 1 I Ll The 

WHICH IS, J £ff - 2 [ 2{ y a _ 2ta) 1- 2 (y a +2t a ) ^ V a 2(V a '-2i a ) 2(V a '+2t a ) V Q ' J ' 1 Ile 

exchange coupling here is antiferromagnetic and is estimated to be (j^J t jr r ' exch 'Jt 2 ^) = 
2.8. It vanishes for U = oo as well as for t a = 0. 

(iv) < r z m >=< r z n >= —1/2 and < t^t* >= 1/4, where the averages are taken over the 
state L\ na l) n JO > (Fig. |3| (e)). Here, the effective exchange coupling turns out to be 
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the same as that of (iii). 

Case VI: A Phase with one Ladder Completely Charge Ordered and the Nearest Neigh- 
bour Disordered: In this case, we have two possibilities, again depending on the state 
through which one computes the averages over the charge part in the effective hamilto- 
nian. In both the cases, one has, < >= =< T y m >, < >= =< > and 
< t^t^ >= =< r^j-y >=< r^jn >. Besides, one has, 

(i) < r z m >= 1/2, < >= —1/2 and < t^t* >= —1/4, where the averages are 

taken over the state a* n JO > (Fig. |3] (f)). All the effective exchange cou- 
plings in this case turn out to be non-zero and hence j™ t j; r > exch j s obtained as, 

jinter,exch _ -t\ h r i l \ j_ \ \ \ l_n 

" 2 l2(V a +2t a ) ~T 2V a (U+2t a ) U (U'+2t a ) U> 2{V^+2t a ) 2V a 'J- 

(ii) < r? n >= —1/2, < >= —1/2 and < t^t* >= 1/4 where the averages are due to 

the state L\ na (v n Q |0 > (Fig. |3| (g)). The effective exchange coupling which is due to 

TT inter(t) -i TT inter(s3) /,i , .1 1 , TT inter(sl) ^ TT inter(s2) . , \ . 

n e ff and H e ff (the contributions due to -n e // and -n e // vanish) is given 

1 jinter,exch _ -t\ h r 1 1 1 ]_-, 

y > e ff 2 \-2(Va+2t a ) ^ 2V a 2(V a '+2t a ) 2V a 'J- 

In both of the above mentioned cases, the exchange coupling becomes antiferromagnetic 
and is estimated to be ^j^f r > exch l /t 2 ab ) = 2.26 for case (i) and 1.65 for (ii). However, it 
vanishes in both cases for U = 00 irrespective of whether V a ^> 2t a or V a <C 2t a . 

IV. DISCUSSION AND CONCLUSION 

We have derived the effective spin-charge hamiltonian for a' -N 3^/2^5 for both intra- 
ladder and inter-ladder exchange. We find a rich structure as a function of possible realiza- 
tion of the microscopic parameters. We have found several, in part unexpected, results. 

(i) The effective magnetic exchange along the ladder decreases with increasing charge order- 
ing (of zig-zag type). For complete charge ordering, the magnetic exchange becomes 
ferromagnetic and quite large in magnitude. 
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(ii) The effective inter-ladder magnetic exchange in between two given rungs of a charge- 

ordered and a charge-disordered ladder changes (only) by a factor of two when the 
charge- density- wave is shifted by a lattice constant along a (compare Fig. |] (f) and 

(g))- 

(iii) There are novel terms of type t^t* in the effective charge-charge inter-ladder interac- 
tion. These terms, which are however rather small in magnitude, could in principle 
stabilize a mixed charge-order configurations like the one illustrated in Fig. |3] (f) and 

(g)- 

As a consequence of (i) and (ii), the proposed frustrated spin-cluster model by de Boer et 
a/.El seems to be unlikely, since (a) no indications of ferromagnetic couplings have been found 
experimentally and (b) the coupling in between the one-dimensional spin-cluster chains of 
RefB should be, as a consequence of (ii), rather strong. This conclusion is consistent with a 
recent study of the frustrated spin-cluster model by DMRG and a cluster-operator theory0. 

Let us note that the change in sign of the effective intra-ladder magnetic exhange shown 
in Fig. [| cannot be described accurately by a perturbation expansion in Near to the 
singularity at V a — 2t a the perturbation expansion breaks down and the effective intra- 
ladder spin-hamiltonian becomes long-ranged. 

Ultimately, the reason for the ferromagnetic intra-ladder coupling for the zig-zag ordering 
found in the above calculation lies in the fact that the charge-ordered state is not the ground- 
state of H . In fact, the charge-ordered state would be stabilized by H' Q in the case of 
a large effective inter-rung Coulomb-coupling Vl — V b " (see Eq. (|13"D). For large values of 
V{ — VI' ^> 2t a , the perturbation expansion would yield (e.g. in mean-field approximation for 
H' ) an antiferromagnetic intra-ladder spin-spin coupling. We did not show this calculation 
here, since our estimated value for the inter-rung Coulomb repulsion Vl — V£ ~ 0.1 eV is 
about one-order of magnitude too small in order to do the job. We therefore believe that 
this small value of Vl — Vl' indicates (a) the importance of elastic effects for the stabilization 
of the observed phase transition at T c = 34 K and (b) that the degree of charge ordering is 
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far from complete. This is consistent with the proposal of about 20% charge orderin; 
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FIGURES 




t V 

t b ,V b ' 



^ab'Xib 



FIG. 1. Schematic structure of a'-NaV^Os where the open circles stands for the Vanadium 
sites. Different parameters are also shown. 






(a) (b) (c) 

FIG. 2. Possible (a) disordered and (b), (c) completely charge ordered (zig-zag) configurations 



in a single ladder. The filled, shaded and open circles denote V +4 , F +4 - 5 and V +5 sites respectively. 
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FIG. 3. Possible (a) disordered; (b), (c), (d), (e) completely charge ordered (zig-zag) and (f), 
(g) ordered-disordered configurations in between two adjacent ladders. 
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(V/2t a ) 

FIG. 4. Variations of the effective intra-ladder exchange coupling Jfff 1 in units of (tl/t a ) 
with respect to the parameter (V a /2t a ) (for U = oo). Note that the perturbation expansion breaks 
down at V a = 2t a . 
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